&INPUT
     LABEL  = 'Two atoms in single trap or adjacent traps with ITYPE=25',
     URED   =   52.547475,
     JTOTL  =   0,   JTOTU  =   0,    IBFIX  =  1,  ! only the block with atom-molecule M=0 is meaningful for trap problem
     RMIN = 0.0,     RMAX = $RMAX,   RMID = 800.0,   RMATCH = $proppar,
     IPROPS =   8,   DRS    =   20.,
     IPROPL =   9,   DRL    =   20.,  TOLHIL = 1.D-6,
     ADIAMN = .FALSE.,
     WKBMN  = .FALSE.,
     BCYOMN = -0.00155038759D0, !-1/a is -1/645
     BCYCMN = -0.00155038759D0,
     !BCYOMN = -1.00D15,         ! close to infinity for harmonic trap without contact interaction
     !BCYCMN = -1.00D15,
     EMIN   = -0.2,  EMAX   =   0.5,  EUNITS = 3,
     DTOL   =  1.D-7,
     IPRINT =   6,
  /

  &BASIS
     ITYPE   = 25,
     JMAX    = 24,    KMAX   = 999,
     JMIN    = 0,
     JSTEP   = 1,       ! JSTEP = 2 for single trap, JSTEP = 1 (default, so can comment) for separated traps
     A(2)    = 1.,      ! Value not used, but non-zero value needed to prevent reading in levels
     ISYM    = 1,  0,   ! Limit basis set to even k, even epsilon * (-1)^k
     BCT     = .TRUE.,  ! Use centrifugal term based on j not JTOT
  /

  &POTL
     RM     =  0.529177210903D0,   EPSIL  = 1.0,
     MXLAM  =  4, LAMBDA =  0,0,            1,0,    2,0,    2,2,
                  NTERM  =  4*-1,
     CFLAG  = 1,        ! Treat potential coefficients as multiplying spherical harmonics C not Y
  /
  &TRAP
   !WRELX=100.D0, WRELY=100.D0, WRELZ=102.D0,
   !WRELX=102.D0, WRELY=100.D0, WRELZ=100.D0,
   WRELX=$omega_x, WRELY=$omega_y, WRELZ=$omega_z,
   DELTAZ = $deltaz,
  /

